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Abstract 

Techniques from effective field theory are applied to nuclear rotation. This approach 
exploits the spontaneous breaking of rotational symmetry and the separation of scale 
between low-energy Nambu-Goldstone rotational modes and high-energy vibrational and 
nucleonic degrees of freedom. A power counting is established and the Hamiltonian is 
constructed at next-to-leading order. 
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1. Introduction 

Collective states are the low-lying excitations of heavy nuclei. One of the hallmarks 
of deformed nuclei are rotational bands of the form 

E(l)&A[l(l+l)-I (I + l)] . (1) 

Here, I > Iq denotes the angular momentum, A is a constant determined by fit to 
data, and Iq is the spin of the band head under consideration. Rotational states are 
the lowest-lying excitations in open-shell heavy nuclei. In rotational even-even nuclei in 
the rare-earth region, for instance, the I = 2 state of the (Iq = 0) ground-state band 
has an excitation energy of about 80 keV while the vibrational Iq — 2 band head is at 
about 1 MeV of excitation energy. Our understanding of nuclear rotation is based on 
the ground-breaking papers by Bohr and by Bohr and Mottelson Here, collective 
nuclear excitations are modeled in terms of quadrupole vibrations of the nuclear surface. 
Rotational nuclei are intrinsically deformed, i.e. they exhibit a mostly axially symmetric 
static deformation in the co-rotating intrinsic reference frame, giving rise to rotational 
bands of the form of Eq. ([1]). The Bohr-Mottelson model has been extended to include 
nuclei with tri-axial deformation within the asymmetric rotor model Rotational nu- 
clear states are also described within the variable moment of inertia model Q, and the 
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general collective model 0, 0, @] ■ The interacting boson model ^S] , an algebraic model 
based on s-wave and <i-wave bosons, is also widely employed for the description of col- 
lective nuclear phenomena. While these models are very appealing due to their physical 
motivation and due to their mathematical beauty, it is difficult to systematically extend 
them, to gauge their limitations, or to compute results with reliable error estimates. Fur- 
thermore, it is non-trivial to keep generalizations of collective models computationally 
tractable This is due to the difficulties posed by the linear realization of rotational 

symmetry. 

The present paper attempts to overcome these limitations. It proposes a model- 
independent description of rotational nuclei that is based on an effective theory (EFT) . It 
treats even-even, odd mass and odd-odd nuclei on an equal footing. In recent years, EFTs 
enjoyed considerable popularity and success in low-energy nuclear structure. Examples 
are the a ppli cation of pion-less EFT to few-body systems [Tljj , dilute Fermi gases with 
repulsive Il2| and attractive [l3l [l4. 15 1 interactions, nuclear interactions based on chiral 
EFT [lj| Tl7l [HI , and the description of halo and cluster states within an EFT [l9|, 
[20l |. Effective field theories exploit a separation of scale, and provide us with a model 
independent description of physical phenomena based on observables 2l|, [22j . They often 
exhibit an impressive efficiency as highlighted by analytical results and economical means 
of calculations. 

Within an effective theory, one proceeds as follows. First, the relevant low-energy 
degrees of freedom have to be identified. For even-even nuclei in the rare earth re- 
gion, for instance, the spins and parities of low-lying states can be explained in terms of 
quadrupole degrees of freedom. Second, the relevant symmetries (and pattern of sponta- 
neous symmetry breaking) have to be identified. For atomic nuclei, the Hamiltonian is 
invariant under rotations. Furthermore, the concept of an intrinsically deformed ground 
state corresponds to the spontaneous breaking of rotational symmetry. (See Bohr's Nobel 
lecture [23j], or Mottelson's lectures in Les Houches [13] ■) Indeed, the rotational spec- 
trum |T|) consists of the low-lying excitations associated with the spontaneous breakdown 
of rotational symmetry. In an infinite system, the corresponding excitations are mass- 
less Nambu-Goldstone bosons. In a finite system such as the atomic nucleus, there is - 
strictly speaking - no spontaneous symmetry breaking and thus no Nambu-Goldstone 
boson. However, as we will see below, the Nambu-Goldstone modes generate the low- 
energetic discrete excitations ((T|) upon the quantization of the finite system 25j. Third, 
we have to identify (and exploit) a separation of scales and introduce a power counting. 
In the case of deformed even-even nuclei in the rare earth region, the separation of scales 
is the separation between, e.g. the rotational J 71 ' = 2 + state (at several tens of keV of en- 
ergy) and the low- lying vibrational 2 + state (at about 1 MeV of excitation energy). Below 
the vibrational threshold, the physics can be described purely within Nambu-Goldstone 
modes. Above the threshold, the vibrations have explicitly to be taken into account, 
and one might also consider the coupling to even higher energetic nucleonic degrees of 
freedom (with a corresponding energy scale of a few to several MeV). For even-even 
nuclei, the phenomenological models jj, 0, @, S H| practically include the physics of the 
Nambu-Goldstone modes and the vibrational degrees of freedom. However, the models 
do not further exploit the separation of scales, and they do not present a power counting 
that would allow for a systematic extension. 

We want to apply the tools of effective field theory to nuclear rotation. The concept of 
spontaneous symmetry breaking is central to this approach. The ground breaking papers 
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by Weinberg |26|, by Coleman, Wess and Zumino [27j . and by Callan, Coleman, Wess 
and Zumino [28[ describe the construction of low-energy Lagrangians in the presence 
of spontaneous symmetry breaking. Leutwyler applied and extended this approach to 
non-relativistic Lagrangian s [29j j , and excellent description of the approach can be found 
in Refs. [il EH, 
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33| . As we will see below, the Nambu-Goldstone fields employed 
in the spontaneous breakdown of rotational symmetry only depend on time (and not on 
space). This is due the finite size of the atomic nucleus. Thus, we do not deal with a field 
theory but rather with quantum mechanics. For this reason, this approach is an effective 
theory (and not an EFT) . Note that nuclear rotation and the spontaneous breakdown of 
rotational symmetry has been addressed within a field theoretical approach by Fujikawa 
and Ui [34] . These authors succeeded in linking nuclear rotation to the Higgs mechanism, 
but they did not pursue the systematic construction of low-energy Lagrangians. 

There are - of course - microscopic approaches to nuclear rotation (see the recent re- 
35| and references therein), and the microscopic computation of the parameters of 
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collective models is a long-standing 
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and interesting problem [37|, Sialic 
of the present paper is not to extract a collective model from an underlying microscopic 
Hamiltonian, but rather to construct a low-energy collective Hamiltonian within an effec- 
tive theory. Apart from symmetry principles, no details of the microscopic Hamiltonian 
are needed for this task. 

This paper is organized as follows. Section [5] presents the derivation of the low-energy 
Lagrangians and Hamiltonians that govern the dynamics of axially deformed nuclei, i.e. 
the spontaneous breakdown of SO(3) to SO(2) is the central concept. In Sect. [3] we 
consider the physics of Nambu-Goldstone modes that result from the symmetry breaking. 
In Sect. S]we study the coupling of vibrations to the Nambu-Goldstone modes. Fermions 
and Nambu-Goldstone modes are coupled in Sect. [SJ The main results of this paper are 
summarized in Sect. [HI Some technical details are presented in the Appendix. 



2. Construction of low-energy Lagrangians 

Let us consider a system with a continuous symmetry such as invariance under ro- 
tations. In this case, the ground state has definite angular momentum. Spontaneous 
symmetry breaking happens when an arbitrarily small symmetry-breaking perturbation 
yields a deformed ground state, i.e. a state with no definite angular momentum (see, 
e.g., Ref. (HI). Clearly, spontaneous symmetry breaking can only take place in infinite 
systems as only these can exhibit gap-less excitations. 

Ferromagnets and antiferromagnets are well known examples for the spontaneous 
breaking of rotational symmetry. In these system, the ground state exhibits a finite 
magnetization and a long-range order of staggered magnetization, respectively. Thus, 
the ground state does not exhibit the full rotational symmetry, but is only invariant 
under rotations around the axis of magnetization. This is the spontaneous breaking of 
rotational symmetry down to axial symmetry. In the case of the ferromagnet, ground 
states with different orientations of the magnetization are inequivalcnt. The effective 
field theory for magnets has been derived by Leutwyler j29|. 

In systems with a finite number of degrees of freedom, such as atomic nuclei, spon- 
taneous symmetry breaking does not occur in a strict sense since arbitrarily small per- 
turbations do not lead to a deformation of the ground state. However, for nuclei in the 
rare earth region, perturbations of the size of a few tens of keV can mix states with 
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different angular momenta. This is a particularly low energy scale compared to other 
excitations, and it is two orders of magnitude smaller than the nucleon separation en- 
ergy. In a technical sense, spontaneous symmetry breaking does not take place because 
differently oriented ground states of intrinsically deformed nuclei are unitary equivalent 
to each other. A superposition (i.e. the projection) of these states creates states with 
good angular momentum that are invariant under the full rotation group |25| . This 
phenomenon is well established in mean field calculations (see, e.g., Refs. |42l. I43L 144]]). 
However, the small energy scale necessary to induce a symmetry breaking justifies and 
motivates us to apply the ideas of spontaneous symmetry breaking to this case. 

For axially deformed nuclei, the rotational symmetry is spontaneously broken. The 
ground state ip^ is only invariant under operations of the subgroup % = 5*0(2), i.e. 

/i^ (0) = ip m for h £ H , (2) 

but it is not invariant under general operations of the full rotation group Q = 50(3) of 
the Hamiltonian. Any rotation g £ Q can be decomposed as g — gh with g £ Q and 
h £ %. Two rotations g = gh and g' = gh' with h,h' £ H that differ from each other by 
an element of the subgroup H yield the same state when applied to the ground state ip^ . 
Thus, such group elements must be identified, and they form an equivalence class. This 
equivalence class is the coset Q/H. The coset SO(3)/SO(2) can thus be used to describe 
the low-energy degrees of freedom which change the orientation of the ground state if>^°\ 
The corresponding degrees of freedom are Nambu-Goldstone modes. In infinite systems, 
these modes have zero mass and thus zero energy in the limit of vanishing momenta. We 
will see below that the Nambu-Goldstone modes of deformed nuclei generate rotational 
bands upon quantization. 

The Nambu-Goldstone modes parameterize the coset 50(3) / 50(2), and we need to 
work out the basic expressions from which rotationally invariant Lagrangians can be 
constructed. This problem was solved for the spontaneous breaking of chiral symmetry 
by Weinberg [26| , and for general Lie groups by Coleman, Wess and Zumino (23] , and by 
Callan, Coleman, Wess and Zumino [28|. For detailed reviews of this matter, the reader 
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Let us parameterize a rotation r £ 50(3) by the three Euler angles a, ft, and 7 as [45[ 

r(a,/3, 7 ) = e -<aJ. e -i/JJ Se -<7-?. . (3) 

Here, k = x,y,z denote the appropriate components of the angular momentum 
operator (i.e. the generators of the Lie group 50(3)). These operators fulfill the usual 
commutation relations 

[Jj, Jk] ^i^SjkiJi ■ (4) 
I 

Let us assume that the subgroup TL = 50(2) consists of the rotations ^1(7) = cxp (— 17 J z ), 
i.e. its generator is J z . This implies that the ground state ip^ has a finite expectation 
value in the z-direction. Thus, any two rotations r(a, ft, 7) that differ by the Euler angle 
7 from each other yield the same state when acting on the ground state. The coset 
50(3)/50(2) thus consists of the rotations 

g(a,ft) = e- iaJ *e- i ? J y , (5) 



and the Euler angles a and (3 are the degrees of freedom of the Nambu-Goldstone modes. 

The dynamics of the Nambu-Goldstone modes is determined by the time derivative 
dtg{pt,P). It is simpler to compute 

,g- 1 (a,/3)9 t5 (a,/3) = iE x J x + iE y J y + iE z J z , (6) 

as this is an element of the Lie algebra of the group Q — SO (3). This defines the functions 
E k , k = x,y,z. 

We want to construct Lagrangians that are invariant under rotations and therefore 
need to study the transformation properties of the functions E x , E y , and E z . Let us 
act with a rotation r £ Q onto a rotation g £ Q /%. This yields rg which again can be 
decomposed into a product of two rotations g £ Q /T-L and h £ H 

rg = g{g,r)h(g,r) . (7) 

Solving for g yields 

9(9, r) = rgh- 1 {g 1 r) . (8) 
The rotations g and h are nonlinear (and complicated) functions of the three Euler angles 
that define r and the two Nambu-Goldstone modes that define g. We have 

g^dtg = (rg^dtirg) 

= h- l (gd t y)h + h- l d t h . (9) 

In this derivation we used that r is time independent, and we employed Eq. ([7]). We 
solve for (gd t g) and find 

g-'dtg = Kg^dtg)^ 1 - (d t h)h^ . (10) 

Similar to Eq. ([5]) we again have 

g~ 1 d t g = iE x j x + iE y J y + iE z J z , (11) 

as this is an element of the Lie algebra of 50(3). We employ Eq. (jlll) and Eq. ^ on the 
left and right-hand side of Eq. (ITU1) . respectively, and observe that hJx^yh^ 1 is a linear 
combination of J x . y while /iJ 2 ft, _1 = J z . Note that the term (dtfyh -1 on the right-hand 
side of Eq. (flO|) is also proportional to J z . Thus, 

E X J X + EyJy — E x hJ x h 1 + EyhJyh 1 , (12) 
E Z J Z = E z J z -i{d t h)h^ . (13) 

These equations show that under a general rotation r, the functions E x and E y transform 
as the x and y components of a vector under the rotation h around the z-axis, while E z 
transforms as a gauge field. Indeed, 

h(d t - lEj^h- 1 =d t - iE z j z . (14) 

Let us further illuminate these derivations. We express the rotation h as h = 
exp (— i"f(g, r) J z ). Here the angle 7(3, r) is a (complicated) function of the Nambu- 
Goldstone modes a and (3 of the rotation g and the three angles that parameterize the 
rotation r. We employ the transformation (|12|) and find 



E x \ = ( cos 7 sin 7 \ / E x 
E y ) \ -sin 7 cos 7 J \ E y 
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(15) 



while Eq. (JT3J) yields 

E z = E z - 7 . (16) 

Note that Eq. (fT5j) implies that 

E ±1 =E x t iE v (17) 

transforms as the upper and lower components of a vector (a spherical tensor of degree 
one), i.e. 

e llh E ±1 = e ±n E ±1 . (18) 

A Lagrangian consisting of combinations of the functions E x , and E y that are formally 
invariant under the subgroup H = 50(2) will thus be invariant under general rotations 
r 6 Q = 50(3) [30|. This is all we need for the construction of Lagrangians involving 
the Nambu-Goldstone modes. Equation (flTj)) shows that the function E z is not invariant 
under a rotation, but merely changes by a total derivative. The function E z will be 
important in cases where time reversal invariance is broken by the ground state (i.e. for 
nuclei with a finite ground-state spin). 

Let us also consider the presence of a field tp that describes physics at a higher energy 
scale. The key is [30j to rewrite 

il> = g(a,P)4,. (19) 

Here, g is the element ([5]) of the coset Q/T-L = SO(3)/SO(2), and (p is defined in terms 
of ip and g. Due to the particular choice (fT9"|). a rotation r € Q transforms ip — > ip as 

ip = rip = rgcp = gh(p ■ (20) 

Here, we used Eq. (0. By definition (compare to Eq. (fTU)) ). we must also have ip = g<f>. 

The comparison with Eq. (|2"0|) shows that under a rotation, the field <p transforms as 
cp —> <p with 

4> = h4> . (21) 

The time derivative of the field 

d t 4>=(d t h)<j> + hd t <t> (22) 
thus transforms as a gauge field. Employing Eq. (|13[) we have 

h(dt-iEj z S }<f>=(d t -iEj z )t- (23) 
Thus, the covariant derivative 

D t = d t - iE z J z , (24) 

when acting onto <p, transforms properly under rotations and can be employed in the 
construction of rotationally invariant Lagrangians for the Nambu-Goldstone modes and 
the field tp. A Lagrangian consisting of of E x , E y , cf>, and D t (p that is formally invariant 
under rotations of the subgroup T-L = SO(2) will be invariant under the full action of the 
group Q — SO(3). We now recognize the advantage of employing a nonlinear realization 
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of the rotational symmetry. While the derivation of the basic building blocks E x , E y , <f>, 
and D t (j) and their transformation properties are somewhat more complicated than for the 
well-known linear representations, the construction of rotationally invariant Lagrangians 
can be achieved by constructing Lagrangians that (at first sight) only appear to exhibit 
axial symmetry. This will make the resulting effective theory computationally tractable. 
This is in contrast to algebraic models which are linear representations of the rotational 
symmetry and employ dynamical symmetries to remain computationally feasible [9] . 

All that remains is to actually compute the functions Ek, k — x,y,z. We use the 
Baker-Campbell-Hausdorff formula and compute the left-hand-side of Eq. ([5]) directly 
from the expression ([5]). This yields 

E x — a sin j3 , 

Ey = ~P, 

E z = -dcos/3 . (25) 

From a practical point of view, these expressions are the main result of this Section. 
Note that the components can be viewed as angular velocities. 

The invariance under rotation implies that angular momentum is a conserved quan- 
tity. To become more familiar with the nonlinear realization of rotational symmetry, we 
compute this conserved quantity via Noether's theorem. Under infinitesimal rotations 
around the k axis (k = x, y, z) by an angle Su>k, the Euler angles a and /? change by the 
infinitesimal amounts Sa, and Sj3, respectively. One finds (see |Appcndix C] for details) 



5a 
5p 




(26) 



with 



jQ. _ ( — cot/Scosa —cot/? sin a 1 \ 
\ — sin a cos a J 

The Lagrangian L(a,/3,j3) of the Nambu-Goldstone modes is invariant under rota- 
tions. We apply Noether's theorem and find the conserved quantities (see | Appendix D| 
for details) 

/ CL \ 

= M T ( Pa I . (28) 




Here, 



Pa 



p * - £' (29) 

to = § (30) 
op 
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are the momenta conjugate to a and /3, respectively. Thus, we have 



Q x — —p/3 sin a — p a cot j3 cos a , 
Qy = P/3 cos a — p a cot (3 sin a , 

<3z = Pq , (31) 
and the total angular momentum squared is 

2 

Q 2 ^Ql + Q 2 y + Q 2 z =p} + ^f- . (32) 

sin /3 

Clearly, Q 2 is the squared angular momentum of a rotor, and we will see in Sect. [3] that 
the leading order Hamiltonian of the Nambu-Goldstone modes is proportional to this 
quantity. 

In what follows, we will continue to employ the Euler angles a and /3 as the Nambu- 
Goldstone modes. This is a convenient but arbitrary choice. The algebraic transfor- 
mation laws derived in this subsection are independent of this particular choice as it 
only depends on the pattern of the symmetry breaking SO(3) —> 50(2). It is only the 
explicit expressions (I25p for the functions E x , E y , and E g , respectively, that depend on 
this parameterization of the coset. 

Let us briefly contrast the effective theory based on nonlinear realizations with the ear- 
lier phenomenological approaches. The Bohr Hamiltonian [l[ and its generalizations [&] 
model the collective states of atomic nuclei by surface vibrations. In leading order, these 
are quadrupole phonons. Within this approach, one transforms from the laboratory 
coordinate system to the co-rotating (or body fixed) coordinate system. Thereby one 
introduces three Euler angles and the two moments of inertia that identify an axially de- 
formed nucleus as the relevant degrees of freedom. Within the effective theory presented 
in this paper, different variables are chosen. In the case of axially deformed nuclei, only 
two Euler angles are relevant for the description of low-energy excitations. As we will see, 
the nonlinear realization naturally corresponds to the spontaneously broken rotational 
symmetry and facilitates the introduction of a power counting. 

3. Nambu-Goldstone modes 

3.1. Even-even nuclei 

Let us assume that the ground state is invariant under time reversal (as is the case 
for even-even nuclei). The simplest Lagrangian is second order in the time derivative, 
and below we will see that this Lagrangian is indeed of leading order. In leading order 
(LO), the Lagrangian of the Nambu-Goldstone modes thus is 

4o } = Y {El + El) = ^ (> + 6? sin 2 /j) . (33) 

Here Co is a low-energy constant to be determined by fit to data. A Legendre transfor- 
mation yields the Hamiltonian 



Here, we employed the canonical momenta pp = dL/d/3 and p a = dL/ da. The Hamilto- 
nian (|34|) is obviously proportional to the squared angular momentum (|32[) . The quan- 
tization in curvilinear coordinates is well known [46}, see |Appendix A| for details. This 
yields the Hamiltonian (we use units where % = 1) 

H = ~f^+oot^ + ^-^) . (35) 
2C V sin p J 

Comparison with the classical Hamiltonian (|3"4")l thus yields 

V% = — =— sin p~dp , 

p Q = -id a . (36) 
The eigenfunctions of the Hamiltonian Q35[) are spherical harmonics, i.e. 

iTY im (/3,a) = l ^±^Y lm {p,a) , (37) 

and I = 0, 1, 2, . . . thus yields a rotational spectrum. Note that the low-energy constant 
Co is, within the collective model, associated with the moment of inertia. 

For even-even nuclei, only even values of I are permitted for the ground-state band. 
This can be understood as follows. The ground state of even-even nuclei is not only 
invariant under rotations h £ H = SO (2) but also under discrete operations such as 
parity, or a rotation about 7r around an axis perpendicular to axial symmetry axis. 
Thus, we must augment the invariant subgroup % by these discrete operations, and this 
modifies the coset accordingly. As a result, the angles (/3,a) and (tt — f3,ir + a) have to 
be identified, and this limits the values of I to even numbers [2] ■ 

3.2. Odd-mass and odd-odd nuclei 

Odd nuclei have a half-integer ground-state spin, while odd-odd nuclei can also exhibit 
a nonzero integer ground-state spin. Thus, the ground state is not invariant under time 
reversal. This modifies the low-energy effective Lagrangian as terms that are not invariant 
under time reversal need to be considered. Such terms consist of only one time derivative, 
and we need to consider the functions (|23|) . The transformation properties (IT21 and (fl~3"|) 
show that the functions E x , E y , and E z are not invariant under rotations. However, the 
function E z only changes by a total derivative, see Eq. (flT)]) . Thus, the action changes by 
an irrelevant phase. In quantum field theory, such a function is known as a Wess-Zumino 
term. In our case, the Wess-Zumino term L\vz is 

L wz = qE z — ~qa cos (3 . (38) 

Recall that a low-energy Lagrangian can be understood as resulting from integrating out 



high-energy fermion modes in a more fundamental Lagrangian 21|. In the case that 
the considered fermion system consists of an odd number of fcrmions, or has a finite 
ground-state spin, the resulting low-energy Lagrangian must reflect this behavior. Thus, 
the appearance of a corresponding symmetry-breaking term is unavoidable. For details, 
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the reader is refered to Refs. [43, Hif- In Sect. [3 we will couple fermions to the Nambu- 
Goldstone modes and see that the relevant term is indeed proportional to E z , as stated 
in Eq. (|55|. 

Let us consider the transformation properties of the Wess-Zumino term. Under a 
rotation by Slo^ around the k axis (k = x,y,z), the Wess-Zumino Lagrangian L\vz 
changes by (see | Appendix D for details) 



and this is obviously a total derivative. The application of Noether's theorem yields the 
conserved quantities 

COS Oi 

Qx = ; — — o — pg sin a — p a cot B cos a. , 

sin p 

„ sin a 

Q v = ; — -q + pa cos a — p a cot p sm a , 

smp 

Qz = Pa , (40) 

which are the components of the angular momentum. The total angular momentum 
squared is 

Q 2 = Q 2 X + Q 2 y + Ql = pl + -^-{p a + q cos pf + q 2 . (41) 

sin p 

With the Wess-Zumino term added, the leading order Lagrangian becomes 



Llo — L LO + Lwz 



2 



(j3 2 + a 2 sin 2 p\ - qa cos 8 . (42) 



The Legendre transformation yields the classical Hamiltonian 

P 2 ( Pa + g cos B) 2 

= ^7=r H , 2 • (43) 

2Go 2Go sm p 

The comparison with the angular momentum (|4"TT) shows that -Hlo = {Q 2 — ? 2 )/(2Co). 
We employ the quantization Q36p. and obtain 



1 

#lo = - or ■ a dp sin Bd R + — (-td a + q cos /3) . (44) 

zGq sm p 2Gq sm p 



The eigenfunctions of this Hamiltonian are products e~ iam d l mq (/3). Here denotes the 
"little" Wigner d function |45|, i.e. the Wigner D function is defined as D l mq (a, /3, 7) = 
e~ ma d mq (f3)e~~ i<rt . Details are presented in the |Appcndix E| Thus, 

H w d l mq ((3)e- iam = E ho (q,l)d l mq (B)e-* am , (45) 
and the eigenvalues are 

1(1 + I) -q 2 



/• 1.0 ':-/•/;• . >r , (i 

kl , i(i + i)-M(M + i) (4()i 



2C*o 2Co 
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Here, q ^ must be integer or half integer, while I formally assumes the values I = 
\q\, \q\ + 1, \q\ + 2, . . ., and m = —I, + (For q — 0, I only assumes even values.) 

Apart from the irrelevant constant \q\/ (2Co), the comparison of Eq. (|4l))) with Eq. (JXJ) 
shows that we have to identify the low-energy constant q with the ground-state spin Iq of 
the nucleus under consideration, i.e. \q\ = Iq. Thus, the spin Iq of the ground-state and 
the spacing between the lowest two states in the rotational spectrum fix the low-energy 
constants q and Cq. This unified description of odd and even nuclei within an effective 
theory is very encouraging. The effective Hamiltonian (|4"4")l yields the correct low-energy 
description of axially deformed nuclei. The case of even-even nuclei (and odd-odd nuclei 
with zero ground-state spin) is particularly simple as q = Iq = 0, while odd-mass nuclei or 
odd-odd nuclei with nonzero ground-state spin Iq have a more complicated Hamiltonian 
due to q = Iq > 0. 

3.3. Next-to-leading order 

We need to establish a power counting. Let us denote the low-energy scale associated 
with the Nambu-Goldstone modes as £. Thus, the leading-order energy (|4T))) scales as 
E-lo ~ £i and we get the following estimates 

c ~ r 1 , 

Pp ~p a ~q ~ ^° , 
~ a ~ E XtVt , ~ £■ (47) 

The estimate in the second line of Eq (|47[) is due to the quantization of angular momentum 
(recall that % = 1), and the estimate in the third line of Eq. ([4"T]) follows from the usual 
relationship between velocities and momenta. It is clear that the Wess-Zumino term 
scales as £wz ~ £, and is indeed of leading order. 

The Nambu-Goldstone modes of energy ~ £ are well separated from the high en- 
ergy (breakdown) scale O ^> £ which is associated with degrees of freedom (vibrational, 
pairing or single-particle degrees of freedom) we have omitted from our theory. In an 
effective theory, these omitted terms manifest themselves through interactions between 
the Nambu-Goldstone modes. At next-to-leading order (NLO), higher derivatives of the 
Nambu-Goldstone modes appear and one new term enters the Lagrangian 



£nlo = L LO + ^(E 2 x +E 2 y ) 2 . (48) 

Here, C2 is the corresponding low-energy constant. The ratio C2/C0 has units of 
energy -2 , and is assumed to be due to omitted physics at the breakdown scale. Thus, 
the dimensional analysis yields 



C2 
Cq 



fr 2 , (49) 



and 



2 



7?(£S + *?)~(h1 «1- (50) 



The last equation indicates that the next-to-leading-order correction is suppressed by 
two powers of £/fi compared to the leading-order terms. 
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For the computation of the Hamiltonian at next-to-leading order, we employ the 
conjugate momenta 



<9Lnlo 



(51) 



Pf3 = 



dp 



Pa = 



da 



(C + C 2 (E 2 X +E 2 y )) 



a sin 2 (3 — q cos f3 , 



(52) 



and apply Eq. (|50l) consistently in the Legendre transform. This yields the Hamiltonian 
at next-to-leading order 



Here, Hlo is the leading-order Hamiltonian (l4"3"l) . Thus, the spectrum of the Nambu- 
Goldstone modes becomes 



where -Elo is given in Eq. (|46[) . Thus, the spectrum at next-to- leading order is a poly- 
nomial of degree two in the right-hand side of Eq. ((T|). This is exactly as given by Bohr 
and Mottelson Q ■ The low-energy constant Ci can be determined by fitting the spacing 
between the rotational state with I = Iq + 4 and I = Iq- For deformed rare earth nuclei, 
the ratio £/fi « 10 _1 , and this explains the high quality of many rotational bands, i.e. 
the next-to-leading-order contribution is a small correction for angular momenta I with 
I <C Higher corrections at next-to-next-to- leading order include terms of the form 



and the dimensional analysis yields C4/C2 ~ il~ 2 . Thus, another factor (£/0) 2 « 10 -2 
separates these contributions from the contributions at next-to-leading order. However 
for rotational states with high angular momentum / « 0(f2/£), the description in terms 
of the Nambu-Goldstone modes must break down, as higher order terms become large. 
In this case, one needs to include degrees of freedom hat are higher in energy such as 
vibrations, pairing effects, and nucleonic excitations. This is the subject of the next two 
sections. 

Note finally that the results presented in this section can also be obtained by simpler 
means. An alternative derivation is given in |Appcndix B| 

4. Coupling of quadrupole bosons to the Nambu-Goldstone modes 

In this Section, we couple higher-energetic phonon degrees of freedom to the Nambu- 
Goldstone modes. The appropriate phonons are quadrupole vibrations, as evident from 
the spins and parities of low-energy states. This results from the observation that spectra 
of even-even nuclei exhibit Iq = 2 band heads at the energy scale O. In this section we 
couple the quadrupole phonons (whose ground state spontaneously breaks the rotational 
symmetry) to the Nambu-Goldstone modes. 




(53) 




(54) 




(55) 
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4-1. Quadrupole phonons 

We consider a quadrupole field ip (spin- two boson) with complex components ip^, 
fi = —2, ... ,2. Invariance under time reversal implies 

i>% = (-1)"^ . (56) 

such that we deal with five degrees of freedom. 

As discussed in Section [21 we parameterize the phonon as tp = g<p with g £ Q/% = 
SO(2>) / SO{2). By assumption, the ground state 0(°) spontaneously breaks rotational 
symmetry, but remains invariant under the operations of the subgroup H = SO (2) of 
the rotation group SO(3). Again, we choose J z as the generator of the subgroup %. 
Consequently, the ground-state expectation value of the field fulfills 

<0 (O) ) M = <, (57) 

with j; / 0, and is obviously invariant under rotations around the (arbitrarily chosen) z 
axis. 

We have to account for the fact that the Nambu-Goldstone modes result from the 
symmetry breaking of the quadrupole phonon 0. Thus, we parameterize the five com- 
ponents of ip = gcj) in terms of the two Nambu-Goldstone modes a and j3, and three 
non-Nambu-Goldstone modes as 

/ 02 \ 



0o ■ (58) 


V ^ / 

The three non-Nambu-Goldstone modes are parameterized by a complex component 
02 and the real component </V The choice (|58"]) for the non-Nambu-Goldstone mode 
is appropriate [30( as is "independent" of the Nambu-Goldstone modes, i.e. it is 
orthogonal to an infinitesimal rotation of the ground state induced by the generators Jk, 
k — x,y that do not belong to the Lie algebra of the subgroup H = SO{2) 

t (J fe (0 (o) )) =0 withfc = x,y. (59) 

The "building blocks" for rotationally invariant Lagrangians are D t tp, Dt<f>2> D t (f>-2, 
E\, E-i, ipo, 02, and 0_2 = 4>2- ^ n y Lagrangian in these quantities that is formally in- 
variant under SO(2) will indeed exhibit full rotational invariance because of the nonlinear 
realization of the symmetry. 

4-2. Transformation properties 

Recall the transformation properties and the conserved quantities of our effective 
theory. Under a rotation r, the field transforms as 

<t)^h{ 1 )4> = e- ll - h (f) . (60) 
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Here 7 = 7(0, /3,r) is a complicated angle of the rotation r and the Nambu-Goldstone 
modes. Note that the component 4>$ is invariant under rotations. For infinitesimal 
rotations by angles 5cok around the fc-axis, we have (see |Appcndix C| for details) 



cosa c sma. 
7 = 3-ow s + -— dujy . (61) 
sin /3 sin /3 



)j, = —2— — -Z2 + P^ cos a — p Q cot p sin a 



Note that 7 = for a rotation around the z-axis. We formally apply Noether's theorem 
to the Lagrangian L(a, P, P, <j>, </>) and find the conserved quantities 

cos a 

Q x = —2 —I2 — Pa sin a — p a cot a cos a , 

sm /3 

^ sin a 
sin/3 

Q 2 = Pa ■ (62) 

Here, we decomposed 4>2 = $2r + i<t>2i into its real and imaginary part, respectively, 
denoted the momentum conjugate to 4>2r {4>2i) by pir (P2i)> and employed the angular 
momentum 

h = {4>2rP2i - 4>2iP2r) ■ (63) 

The total angular momentum squared is 

Q 2 = p i p + — s - (pi + 4p a / 2 cos/3 + 4Z|) 
sm p 

= P0 + "T^U (Pi + 2h cos/3) 2 + (2/ 2 ) 2 . (64) 
H sm p 

4-3. Power counting 

Let us consider the power counting. There are two possibly distinct breakdown scales 
for our effective theory. First, the restoration of rotational symmetry at high excitation 
energies signals the breakdown of the effective theory. This scale results from the scales 
£ and f2 that describe the quadrupole phonon ip = gcj>. Second, additional degrees of 
freedom enter at an energy scale A 3> fl, and their effect also needs to be considered. 
(We neglect, for instance, collective phonons of higher multipolarity, pairing, and single- 
particle degrees of freedom.) We will focus here primarily on the restoration of rotational 
symmetry. The effects from not included high-lying degrees of freedom will be neglected 
by formally sending A — > oo. 

The field (|58p spontaneously breaks the rotational symmetry, and we separate the 
vacuum expectation value v from the fluctuating contribution tpo as 

tpo = (f) - v . (65) 

We have the following scaling relations 

u~0o ~ r 1/2 , 

Dm ~ D t <h ~ W 2 , 

= 0o ~ 4>2 ~ V 1 ' 2 , (66) 
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in addition to the relations (|47[) . The finite vacuum expectation value v must clearly be 
associated with the low-energy scale £, while the fluctuations of the field <\> are associ- 
ated with the high-energy scale f2, and are necessarily much smaller than the vacuum 
expectation value. Indeed, we have ipo/v ~<f> 2 /v~ (C/^) 1/2 < 1. 

Let us briefly discuss the case of a finite breakdown scale A ^> ft. Consider as 
an example the term C(j)\. Here, C has dimensions of energy -1 , and it must scale as 
C ~ A -1 . Thus, the contribution C§\ ~ il 2 /A is of next-to-leading order, i.e. it corrects 
the leading-order term that scale as f2. By sending A — > oo, we neglect such terms. 
However, it is clear that a full-fledged effective theory needs to deal with them at some 
point in the power counting. The importance of such terms depends on how the ratio 
that governs the validity of the spontaneous symmetry breaking compares to the 
ratio O/A that governs the relevance of neglected degrees of freedom. 

The construction of low-energy Lagrangians in the presence of a spontaneously broken 
symmetry has been discussed in Section^ with a focus on kinetic terms. We also need to 
discuss the form of admissible potentials V(<f>). A rotationally invariant potential V(<f>) 
consists of expressions involving <j> that are formally invariant under rotations around the 
z axis. Furthermore, the potential must exhibit spontaneous symmetry breaking. We 
expand V — Vlo + Vnlo + In leading order we have 



The low energy constants must scale as ujq ~ 0J2 ~ Q to yield the ground-state expectation 
value (Vlo) ~ SI- 

For the construction of next-to-leading order potential terms, we need to determine 
the breakdown scale for our effective theory. With increasing energy, the fluctuations 
ipo grow in size, and the effective theory breaks down once ipo ss v, since this implies a 
restoration of the spherical symmetry. Likewise, large excitation energies correspond to 
a large amplitude \4>2\ ~ v, again resulting in a restoration of the spherical symmetry 
and in the breakdown of the effective theory. The minimum kinetic energy of a large- 
amplitude field <j) is v~ 2 ~ £, while its potential energy is (Vlo) ~ ^ 2 /f 3> fi< We will 
see that the power counting has to employ the kinetic energy scale. 

Let us make a polynomial expansion of the potential 



Here, k,l are integers, and it is understood that only terms with k + 21 > 2 are being 
summed over. (The leading order terms are k + 21 = 2.) Clearly, the potential exhibits 
spontaneous symmetry breaking and is formally invariant under SO(2). Due to the 
nonlinear realization of the SO(3) symmetry, it is also invariant under SO(3). 

Note that Vkiifio\4>2\ 21 /Vlo is dimensionless. This implies that Vki/ti 2 has dimension 
of energy' - 1+fc / 2 . For the power counting we have to assume that this energy scale has 
to be identified with the kinetic energy scale. Thus, we assume 



^Lo(0) = ^(0O-,) 2 + ^|0 2 | 2 . 



(67) 




(68) 



fc+2Z>2 



(69) 



and find 




(70) 
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This establishes the power counting for the potential. Clearly, for large amplitudes 
¥>o> 1021 ~ £ -1 / 2 each potential term of the sum ([68]) is of order Q, signaling the break 
down of the effective theory. 

Our power counting is valid for well deformed nuclei that are "rigid" rotors , i.e. 
nuclei for which loq,lo2 3> £■ Some deformed nuclei are not in this category. So-called 
7-soft nuclei, for instance, do not exhibit a separation of scale between excitations within 
a rotational band (energy scale £) and the vibration with frequency For these nuclei, 
the power counting is different. 

4-4- Even-even nuclei 

Let us consider the case where the symmetry-breaking ground state is invariant under 
time reversal. The following kinetic terms with two time derivatives are invariant under 
rotations 



(D t <po) = tpo 
D t faD t fa-2 = 



1,2 

2 +4lia^fafa^E z +A\fa\ 2 E 2 
E'l + El = /3 2 + a 2 sin 2 /3 . (71) 
Note that, in leading order, the covariant derivative is simply the usual time derivative. 

4.4.1. Leading order 

In leading order (0(fl)), the Lagrangian is 

1 . 2 • 2 UJ 2 n UJ 2 



f^-^1021 2 . (72) 



This Lagrangian describes harmonic vibrations of the quadrupole degrees of freedom. 
The Legendre transform yields the Hamiltonian 

H LO = \pI + \ {vlr + A) + + ^ + A) ■ (73) 

Here, fa r and fai are real variables that denote the real and imaginary parts of fa, 
respectively, i. e. 



b 2r + ifa t • (74) 



The momenta are defined as 



dLho <9Llo <9Llo /-n 

P0 = -JTT- , Vlr = —7— , P2i = —7— • (75) 

o<Po dfa r dfai 

One might wonder whether the most general leading-order Lagrangian should not have 
dimensionless low-energy constants in front of its kinetic terms. Note, however, that any 
such constants could be absorbed by a redefinition of the variables ipo and fa, and a 
rescaling of the oscillator frequencies ujq and UJ2- 

The leading order part of the Hamiltonian (|73l) describes an axially symmetric har- 
monic oscillator in three dimensions. The canonical quantization rules po —y —id Vo , 
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p 2r — > —idrj, 2r ., and p 2 i — > ~id$ 2i yields the Hamiltonian operator £/lo- It is of advan- 
tage to seek not the Cartesian eigenstates but rather eigenstates which reflect the axial 
symmetry. To his purpose, we write <j) 2 = <p 2 e n and obtain the leading-order Hamiltonian 




Hlo = -\dl + ^ - T ( -<> h - —it,, + „-j T -j j . (70) 

Here, 

h = -t&, (77) 

is the operator corresponding to the classical expression ([55]). The eigenstates and 
eigenenergies fulfill 

HLo\n n 2 h) = E(n ,n 2 ,l 2 )\n n 2 l 2 ) ■ (78) 



The energies are 



£■(",)• ii- 2 . 1 2 ) =-n j no + \) + y(2?i 2 + + ■ (79) 



Here, hq — 0, 1, 2, . . . denotes the quantum number in the tpo-coordinate (i.e. excitation 
along the symmetry axis), n 2 = 0, 1, 2, . . . denotes the radial quantum number associated 
with the radius ip 2 — {4>\ r + fili) 1 ^ 2 , while l 2 = 0, ±1, ±2, . . . denotes azimuthal quantum 
number, i.e. the quantum number of the operator l 2 . Clearly, the energies (|79[) are of 
order f2. 

4- 4-2. Next-to-leading order 

In next-to-leading order (O (£,)), the Lagrangian is 

,2 



V - {E 2 X + E 2 y ) - 4E z lm (j> 2 <j>*) + i?NLO (80) 



-^NLO = -^LO 

with 

-Rnlo = (Z 00 ipl + Z 02 \(f) 2 \ 2 )(pl + (Z 20 (fio + Z^l^jlDtc, 

£ w k \h\ 21 ■ (81) 

fc+2i=3,4 

The terms of the Lagrangian describe the rotations, and the coupling between rota- 
tions and vibrations, while anharmonic corrections to the vibrations are encoded in the re- 
mainder (|8ip. It is assumed that the "wave function renormalizations" Zqq, Z 2 q, Zq 2 , Z 22 
are of order £. Under this assumption, the remainder -Rnlo merely adds perturbative cor- 
rections to the vibrations. As we are mainly interested in the lowest energetic vibrational 
states (the band heads) and the modifications of the rotational bands due to the vibra- 
tional coupling, we do not consider these anharmonicities. A Legendre transformation 
yields the Hamiltonian 

i?NLO = H LO + ^2 [Pi + 7^ (Pa - 2/ 2 COS Pf ) • (82) 



sin 
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The comparison with the angular momentum (1641) shows that 



-Hnlo — Hlo + 7T 1 ? (Q 2 ~~ ^2) • 



(83) 



Note that the identification of the squared vacuum expectation value v 2 with the moment 
of inertia is arbitrary but convenient. We could also have employed the low-energy 
constant Co and avoided any reference to v. The vacuum expectation value v is not an 
observable, and the most general Lagrangian can be expressed in terms of the variable 
fo instead of (f>o = v + ipo- 

Let us turn to the quantization of the Hamiltonian at next-to-leading order. From 
the simple form of the classical Hamiltonian (|83p we expect a rotational band upon each 
vibrational state, with a spin of the band-head equal to 2Z2- It is instructive to derive this 
anticipated result. For the quantization of the Nambu-Goldstone modes we use Eq. ([3"6")) . 
This yields the Hamiltonian operator 



The eigenfunctions of i?NLO — -^lo are again given in terms of the Wigner d function as 
e~ ma d ml2 (j3). The energies at order £ thus are 



Here, the angular momentum I is an integer with I > 2 1 Z2 1 and depends on the I2 value 
of the vibrational state with quantum numbers (no,n, 1%) under consideration. The 
spectrum (|85|) is the well known rotational-vibrational spectrum: upon each vibrational 
state, there is a rotational band. Note that the moment of inertia is, at this order of the 
effective theory, equal for each vibrational state. 

The low-energy constants of the effective theory are the vacuum expectation value v 
(which fixes the moment of inertia), and the frequencies ujq and CJ2 (which determine the 
spacing of the ground state to so-called /3 band with quantum numbers (no = l,ns — 
0, £2 = 0) and the so-called 7 band with quantum numbers (no = 0,n2 = 0, Z2 = 
respectively. If we are only interested in these three bands (as is often the case), there is 
no need to determine the low-energy constants of the anharmonic corrections (|8ip . 

4-4-3- Next-to-next-to-leading order 

At next-to-next-to-leading order (N 2 LO) (0(£ 2 /f2)), new terms enter. There will be 
(relatively uninteresting) terms i?N 2 LO that only depend on the vibrational degrees of 
freedom. More interesting are the additional terms that couple rotations and vibrations. 
These are 




(84) 




(85) 



AL N 2 



N 2 LO 



+ D (El+El)yl+F Q (El + El)yl 



+ 




(86) 
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Here, Do, Di, and Di are dimensionless low-energy constants of order one, while F\, F2, 
and F3 are expected to be of order f2~ 2 . Many of these terms are not considered in the 
generalizations 0, HI, @ of the Bohr Hamiltonian. 

Let us briefly discuss the quantization procedure. As is well known, the quantization 
of a classical Hamiltonian is only without ambiguity in flat space, i.e. when the metric 
exhibits no curvature. For a metric with nonzero curvature, one can think of the system 
as a constrained system, i.e. the motion is constrained to a (curved) hypersurface in a 
higher-dimensional space. Kaplan, Maitra, and Heller [49j showed that the quantization 
of such a constrained system depends on the exact nature of the constraints, i.e. on 
details regarding the implementation of the constraining forces which - in the limit of 
infinite forces - confine the motion to the hypersurface. These authors point out that 
the physical situation (i.e. the implementation of the constraints in a limiting process) 
resolves the ambiguity. Our quantization follows this rule. The leading order vibrational 
Hamiltonian is quantized without ambiguity. When the rotations are coupled to the 
vibrations, we know that in the limit of infinite frequencies ujq and ui2 - the physics 
of Nambu-Goldstone bosons as presented in Sect. [3] must result. Thus, we employ the 
previously derived quantization rules. When considering terms beyond next-to-leading 
order, the quantization rules remain unchanged. 



5. Coupling of fermions to Nambu-Goldstone modes 

For deformed odd-mass nuclei, it is not sufficient to simply include terms that are 
first order in the time derivative (such as Wess-Zumino terms) into the Lagrangian. Such 
nuclei oft exhibit several low-lying band heads with rotational bands upon them. In this 
case, fermionic degrees of freedom must be included in the description, and we need to 
coupled nucleon fields N to Nambu-Goldstone modes. As discussed in Section [2j we 
write 

N = g X , (87) 

where g £ TL — 50(2) as given in Eq. ((3]). The most general rotationally invariant La- 
grangian consists of combination of \, D t x, and the functions E x and E y that is formally 
invariant under rotations around the z-axis. In lowest order, we have the Lagrangian [50] 

L = y {El + El) + x f Tx + x f (idt + E Z J Z ) X + (eJ x + E y J y ) x . (88) 

Here, T is the operator of the kinetic energy (or the single-particle energy) for the fermion. 
The Lagrangian ([55)1 reminds us of the particle-rotor model 0, H[ • The choice of the spin 
q of the fermion field is determined by the nucleus under consideration. On could, of 
course, identify x with a two-spinor x(r) (Then, we would also need to integrate over 
space in the Lagrangian.) However, for heavy nuclei, it might be more adequate to 
identify x with the shell-model orbital (n,l,j,T z ) (denoting radial quantum number, 
orbital angular momentum, total spin, and isospin projection, respectively) that is most 
relevant for the open-shell nucleus under consideration. 

It is insightful to consider the limit of a spin-g fermion with components x = {Xq> 0) T . 
Let us assume that Xq is the ground state of the operator T, with T ~ A ^> £. In this 
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case, the low-energy Lagrangian (|88| for the Nambu-Goldstone modes becomes (we use 

l = y( e - + E v) + <* E * + • ( 89 ) 

Apart from the irrelevant constant (x|T|%), this Lagrangian has the same form as the 
Lagrangian (|42p with the Wess-Zumino term. This example shows how the Wess-Zumino 
term arises from high-energy fermionic degrees of freedom. The detailed discussion of 
fermions coupled to Nambu-Goldstone modes is beyond the scope of this paper. However, 
the techniques we developed in this paper seem to be applicable to this case, too. For a 
nucleus with a finite ground-state spin, one would need to couple the Nambu-Goldstone 
modes to quadrupole bosons and single-particle degrees of freedom. 



6. Summary 

This paper developed an effective theory for deformed nuclei. The approach ex- 
ploits the spontaneous breaking of the rotational symmetry and is based on the separa- 
tion of scale between low-energetic Nambu-Goldstone rotational modes, higher energetic 
quadrupole modes, and single-particle degrees of freedom at much higher energies. The 
nonlinear realization of the rotation group is key to the effective theory, and a power 
counting is established. We derived the Lagrangian (and Hamiltonian) for the Nambu- 
Goldstone modes and the for the quadrupole bosons coupled to Nambu-Goldstone modes 
in next-to-leading order. At this order in the power counting, well known results from 
phenomenological models were rederived in a mo del- independent way. More interesting 
phenomena are expected at next-to-next-to-leading order, as the effective theory predicts 
the appearance of terms that are not employed in the phenomenological models. The ef- 
fective theory treats nuclei with finite spins in their ground states (odd-odd and odd-mass 
nuclei) on equal footing to nuclei with zero ground-states spins (even-even nuclei). 

Acknowledgments 

The author gratefully acknowledges discussions with R. J. Furnstahl, W. Nazarewicz, 
A. Schwenk, and H. A. Weidenmiiller. He thanks the Institut fur Kernphysik, Technische 
Universitat Darmstadt, and the GSI Helmholtzzentrum fur Schwerionenforschung for 
their hospitality. This work was supported by the U.S. Department of Energy under 
Grants Nos. DE-FG02-07ER41529 and DE-FG02-96ER40963, and by the Alexander von 
Humboldt-Stiftung. 

Appendix A. Legendre Transformation 

Let the Lagrangian depend on velocities dtx (x = (xi, . . . , xn)) as follows 

L = - {d t xf Gd t x + k T d t x . (A.l) 

Here G denotes the symmetric mass matrix, T denotes the transpose, and k is a constant 
vector. Then, the Legendre transformation 

» = % « A - 2 > 
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yields the Hamiltonian function 



H = \(p-k] T G-Up-k) . (A.3) 



For the quantization of the Nambu-Goldstone modes (a, (3), we follow [46| and write 
the Lagrangian as a quadratic form 

L=\0,a)G ( I ) . (A.4) 

This defines the (2 x 2) matrix of the metric G. The Hamiltonian becomes (let us use 
units where h = 1) 



H = 1 



2v / detG 



(d^d^G^Vtetdf d g p ^ . (A.5) 



For an in-depth discussion of the quantization of constrained systems, the reader is refered 
to Ref. (H. 



Appendix B. Nambu-Goldstone modes revisited 

The derivations presented in Sect. [5] are more formal than necessary if one is only 
interested in the physics of Nambu-Goldstone modes. Here, we follow Leutwyler [29[ for 
a quicker derivation. The dynamics of the Nambu-Goldstone modes is determined by 
the coset SO(3)/SO(2) which is isomorphic to the two-sphere S 2 . Thus, the Nambu- 
Goldstone modes parameterize the two-sphere, and we can choose the parameterization 

(cos a sin /3 \ 
sin a sin (3 . (B.l) 
cos£ / 

Here, a and (3 are time-dependent variables. The simplest low-energy Lagrangian that 
is invariant under time reversal is 

L = ^(a t n).(a t n)^^(/3 2 + d 2 sin 2 /3) . (B.2) 

This is Eq. (|3"3"|) . One can pursue this direction further and also construct the Wess- 
Zumino term. Details are given in Ref. [29j . 



Appendix C. Transformation properties under rotations 

Let 

g{a,P) = e - ia3 *e- ifs3 * (C.l) 
denote an element of the coset SO{2>) / SO{2) and 

r(a,/3, 7 ) = e-toS'e-V^e-* 3 ' (C.2) 
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be a general rotation parameterized in terms of the three Euler angles. Let 



ft( 7 ) = e-' 1 ^ (C.3) 

be a rotation of the subgroup H = 5*0(2). 

The product rg is again a rotation, and we have 

r{a 2 , p 2 , J2)g(a, /3) = r(a,p~, 7) . (C.4) 

Expressions for (a,/3,7) in terms of the other angles are well known [45| 

cot(<S — 02) = cos /3 2 cot(a + 72) + cot (3 — s ™-@ 2 — 

sm(a + 72) 

cos/3 = cos/3 cos/?2 — sin/3sin/? 2 cos(a + 72) , 

cot7 = cos/3cot(a + 72) + cot/3 2 — sin ft — (C.5) 

sm(a + 72) 

Due to the definitions (|C.1[) - (|C3|) . we also have 

r(aj,j) = g(aj)h(j) . (C.6) 

Under a rotation r(a.2, P2, 72) the element g of the coset transforms as 

g = g(a,P)^g = g(aJ) , (C.7) 

and g = g(g, r) depends on the Nambu-Goldstone modes (a, /3) that parameterize g and 
the rotation angles (02,^2, 72) of r. The rotation thus maps the Nambu-Goldstone modes 
(a,/3) into (a,/3), and one needs to employ the transformations (|C5[) to obtain explicit 
results. These transformation laws are complicated, but we do not need them for the 
construction of Lagrangians that are invariant under rotations (as shown in Subsection^ . 
They do, however, enter the derivation of the conserved quantities, i.e. the components 
of the angular momentum. For this purpose, we need to know the transformation laws 
for rotations by infinitesimal angles 5u>k around the k = x,y,z axes. The corresponding 
rotations are r(— 7r/2, 8ui Xl 7r/2), r(0, 8uj y , 0), and r(Suj z , 0, 0) respectively. Employing the 
transformation laws (|C.5|) we find 

5a 

5/3 

with 




= M \ Suj y , (C.8) 



- _ / — cot /? cos a — cot /? sin a 1 \ . . 

Y — sin a cos a y ' 

This is Eq. (J2BJ). 

We can repeat these considerations for the quadrupole phonons. As shown in Sect.[H 
under rotations r(a 2 , /?2 ? 72)5 the quadrupole modes transform as 

-> /i( 7 )<£ = e-^"4> . (CIO) 
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According to the transformation laws (|C5[) . under infinitesimal rotations by 8u> around 
the x, y, and z axis, the angle 7 becomes 7 = Suj^j, 7 = <5wf^|, and 7 = respectively. 
The application of Eq. (|C.10j) shows that the quadrupole fields (with 2 = 4>ir + i<j>n) 
transform as 



(C.ll) 



5<t>2r ^ 


, / Sui x 




\=n[ 5oj y 


$<Po 1 


* \ Slo z 



where 




" = -2Wr£far -2 s ^%hr . (C.12) 



As expected for the nonlinear realization, ipo is invariant under rotations. We can use 
these transformation laws to apply Noether's theorem. 



Appendix D. Application of Noether's theorem 

Let us recall Noether's theorem for our purposes. Consider a Lagrangian L of coordi- 
nates q v and velocities q v , v = 1, . . . , N. Let us also consider a coordinate transformation 
that depends on parameters LJk, k = 1, . . . , K. For infinitesimal small arguments Slu^ the 
coordinates change by 

K 

5q v = ^M vk 8uj k . (D.l) 

k=l 

The coordinate transformation changes the Lagrangian by the amount 

N 

j OLi , 
I An -I- 



«■- £ ■ < D - 2 > 



i/=i 

We employ the equations of motions 



|^ = ^, v = l,...,N, (D.3) 
dq v dq v 



and find 

N 

dqv 



v=l 



Let us consider the case that only one parameter Lj k is varied. Thus, the Lagrangian 
changes by 

5L = 5u) k d t V —M vk • (D.5) 
tH 9 ^ 
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In our case, the transformations are rotations. In the case that the ground state is 
invariant under time reversal, we have SL = for arbitrary Su>k- Thus, 



N dL N 
® k = 2l-Q^ M »k=2_,P»M vk (D.6) 

V— 1 V— 1 

is a conserved quantity (i.e. dtQk = 0). Here, p v is the momentum conjugate to q v . 
We insert the matrix M into Eq. (|D.6I) and obtain the three components of the angular 
momentum (|31l) . 

We can repeat these considerations when quadrupole phonons are coupled to the 
Nambu-Goldstone modes. The application of Noether's theorem yields 

(D.7) 




Thus (h = &2rP2i ~ 4>2iP2r), 

cos a 

Q x = —2 — — -<2 — pa sm a — p a cot p cos a , 
sin p 

sin a 

Q v — —2— — -t2 + vr cos a — p a cot p sma , 
sm p 

Q 2 = Pa ■ (D.8) 

This is Eq. (|52]). 

In the case that the ground state breaks time reversal symmetry, the Lagrangian is 
not invariant under infinitesimal rotations but changes by the amount p9j) due to the 
Wess-Zumino term. Thus, 

SL = {iwz -.(***(=f) + Kft (Sf ) ) ■ CD.») 



This change clearly is a total time derivative, and equating the expressions (ID.5|) and 



(ID.9[) for k = x,y,z yields the conserved quantities (|40 



Appendix E. Solution of differential equation 

We discuss the diagonalization of the Hamiltonians ((33]) and (|84l) . The Nambu- 
Goldstone modes are essentially governed by the Hamiltonian 

11 2 
H = — —-da sin 0da + . , - (-id a + qcos/3) . (E.l) 
sm p sin /3 

For the eigenfunction (a,/3), we make the ansatz 

/i,(a,/3)= 5 ^(/3)e- ima • (E.2) 

Here / is a quantum number that needs to be determined. This yields the eigenvalue 
equation 

j — ^—zdpsinpdp + . \ (m- q cos f3) 2 \ g mq {(3) = E(l,m,q)g mq (/3) . (E.3) 
^ smp sm p J 
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We expand the square, rewrite cos 2 (3 = 1 — sin 2 (3 and find 



-dpsmfidf} + -r^- (m 2 - 2mq cos (3 + q 2 ) \ g l mq (f3) 



sin (3 p sin (3 

= (E(l,m, q )+q 2 )g l mq ((3) . (E.4) 

For the eigenfunctions of this differential operator, we recall that the Wigner D functions 

D l mq (a, f3, 7 ) = e~ ima d l mq ((3)e'^ (E.5) 

(and the functions d l mq ((3) themselves) solve the differential equation (iBj 

(sin/%) + -\- (to 2 - 2m<7Cos/3 + q 2 ) \ D, 1 (a, 
sinp sin p J 

= l{l + l)D l mq {a^ n ) . (E.6) 

Thus, we have to identify the eigenfunctions as q l mq (!3) = d l mq ((3), and the eigenvalue 
is E(l,m,q) = 1(1 + 1) — q 2 . It is well known [45[ that for q — 0, the Wigner function 
D l m0 (a, (3,j) is proportional to the spherical harmonics Yj_ m (/3, a). 
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